Quantum calculations of H 2 — H 2 collisions: from ultracold to thermal energies 
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We present quantum dynamics of collisions between two para-Eb molecules from low (10 -3 K) to 
high collision energies (1 eV). The calculations are carried out using a quantum scattering code that 
solves the time-independent Schrodinger equation in its full dimensionality without any decoupling 
approximations. The six-dimensional potential energy surface for the H4 system developed by 
Boothroyd et al. [J. Chem. Phys. 116, 666 (2002)] is used in the calculations. Elastic, inelastic 
and state-to-state cross sections as well as rate coefficients from T — 1 K to 400 K obtained from 
our calculations are compared with available experimental and theoretical results. Overall, good 
agreement is obtained with previous studies. 

PACS numbers: 



I. INTRODUCTION 

Collisions between molecules are fundamental pro- 
cesses which take place in many areas of physics and 
chemistry. Molecular collisions are prevalent in the in- 
terstellar medium, in the atmospheres of planets, in com- 
bustion chemistry and and in many other chemical and 
industrial processes. Being the simplest neutral molecule 
- molecule system, the H2— H2 system has served as a 
prototype for accurate calculations of tetratomic poten- 
tial energy surfaces (PESs) as well as accurate quan- 
tum dynamics treatment of diatom - diatom collisions. 
Since H2 is the most abundant molecular species in the 
universe, collisions between H2 molecules leading to ex- 
citation and de-excitation of their rovibrational levels 
continue to be a topic of considerable interest in inter- 
stellar chemistry. However, an accurate calculation of 
rovibrational energy transfer rate coefficients in H2— H2 
collisions for temperatures relevant to astrophysical en- 
vironments is still a challenging problem. In recent 
years, due to the success in creating dense samples of ul- 
tracold molecules in magneto-optical traps, atomic and 
molecular collisions in cold and ultracold gases have at- 
tracted considerable experimental and theoretical atten- 
tion 0, i, [3, i, i, @, 0, S|. These unusual systems 
provide a fascinating opportunity to investigate atomic 
and molecular processes at temperatures close to abso- 
lute zero where the collisional outcomes are dramatically 
influenced by quantum effects. Recent studies of atom - 
diatom inelastic and reactive collisions at temperatures 
close to absolute zero have shown that such processes 
may occur with significant rate coefficients at ultracold 
temperatures @, Eft EL HE El HI HE HE HI EE El- 
The H 2 — H 2 system serves as an ideal prototype to inves- 
tigate molecule - molecule collisions at ultracold temper- 
atures. 

Several PESs of the H4 system have been calculated 
in the literature. Some of the earlier studies treated 
the diatomic molecules as a ri gid rotor (RR) as in the 
work of Zarur and Rabitz (ZR) [2fj|, Schaefer and Kohlcr 
(SK) [U, Diep and Johnson (DJ) (HI, and the recent 



study of Patkowski et al. [23||. Some of these PESs have 
been employed in several studies of rotational energy 
transfer in H2— H2 collisions within the r igid-ro tor ap- 



proximation at thermal energie s I2C 
and also at ultralow energies 24 



I26|,l27|,|28|,|2gl 

. Several full- 
dimensional potential surfaces have also been reported 
for the H4 system in recent years. They include the 
PESs of Schwenke [sjj, Aguado, Suarez and Paniagua 
(ASP) [31f. Boothrovd. Martin, Keogh and Peterson 
(BMKP) [32j, and the more recent work of Hinde [33|. 
These PESs have been adopted in a number of quan- 
tum scattering calculations of H?— H2 th at g o beyond 
the rigid rotor approximation [3J, HE HE 1371 [38j . The 
BMKP and Schwenke PESs were recently employed in 
full-dimensional quantum calculations by Pogrebyna et 
al. [3E E3|- They reported the first quantum calculation 
of vibrational relaxation in collisions between H2(i> = 
l,j = 0) and H 2 (w = 0,j = 0) molecules using a 
time-independent quantum mechanical (TIQM) method 
within an angular momentum decoupling approximation 
called the coupled-states approximation (CSA). Lin et 
al. EH used a time-dependent wave packet (TDWP) 
method and the CSA method to study pure rotational 
transitions in H2(u = 0,j = 0) + H2(v = 0,j = 0) 
collisions at thermal energies. Full-dimensional quan- 
tum studies without CSA have also been recently re- 
ported using the TDWP methods El EE El- G&tti 
ct al. [HJ and Otto et al. E3 investigated para-para col- 
lisions involving H-Av = 0,j = 0) + H 2 (u = 0,j = 0), 
and Panda et al. E3 reported ortho-para collisions of 
H 2 (u = 1,0, j = 1) + H 2 (w = 0,1, j = 0). In a re- 
cent Communication [45| , we reported rotational and vi- 
brational relaxation of H2(w = l,j = 0,2) in collisions 
with H.2(u = 0,j = 0,2) at ultralow energies using the 
BMKP surface and a TIQM method that does not in- 
volve any angular momentum decoupling approximation. 
The calculations showed that indistinguishable molecule 
- molecule collisions may involve highly efficient near- 
resonant energy transfer if they are accompanied by si- 
multaneous conservation of the total rotational angular 
momentum and internal energy of the molecules. This 
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mechanism was found to be independent of the initial 
vibrational excitation of the molecules. 

In this article, we present the full quantum dynam- 
ics of rovibrational energy transfer in para-H2 + para- 
H2 system from ultralow to high collision energies and 
compare our results with previous quantum calcula- 
tions l24l [27l. [3 9l. l43j and available experimental re- 
sults [2^ SaEjT^ We employ the BMKP PES and the 
full-dimensional TwoBC - quantum scattering program 
developed by Krems [48| implemented in our previous 
study 45j. The paper is organized as follows: In sec- 



tion II, we give an overview of the TIQM scattering the- 
ory of two 1 E diatomic molecules. In section III, we 
present state-to-state and initial state-selected cross sec- 
tions and rate coefficients for the H 2 (v = l,j = 0) + 
H 2 (u = 0,j = 0) andH 2 ( V = 0,j = 0) + U 2 (v = 0,j = 0) 
collisions. Conclusions are given in section IV. 



T(ri) and T(r 2 ) for each diatomic molecule and the PES 
function, 



V(r u r 2 ,R) = U(r u r 2 ,R) + V(n) + V(f 2 ), 



(2) 



describing the interactions between the four atoms. The 
terms V(fi) and V(f 2 ) in Eq. (J2j) are the interaction en- 
ergy potentials of the two H2 molecules and U(fi,f 2 , R) 
is the interaction energy potential between the two 
molecules, which vanishes at large molecule - molecule 
separations, 



lim U(fi,f 2 ,R) 

R— *oo 



0. 



(3) 



The angular dependence of the interaction potential may 
be expanded as 

U(n,r 2 ,R) = Yl A * ( r i . r 2 . R ) Y * ^ > R ) > ( 4 ) 



II. MOLECULE 



MOLECULE SCATTERING 
THEORY 



The TIQM formalism for scattering of two 1 S di- 
atomic molecules has been described by Takayanagi [4!| , 
Green 15011 . Alexander and DePristo [51(, and Zarur and 
Rabitz [20], using the close-coupling (CC) formalism of 
Arthurs and Dalgarno (52j . A brief account of the the- 
oretical formalism with particular emphasis on distin- 
guishable versus indistinguishable molecule collisions is 
given here. For more details we refer to the original pa- 
pers. 




R 




FIG. 1: Jacobi vectors used to describe molecule 
collisions. 



molecule 



Figure [T] shows the Jacobi vectors employed in the cal- 
culations where f\(r\, r±) and r 2 (r 2 , r 2 ) describe the vec- 
tors joining the two atoms of the two H2 molecules and 
R(R, R) denote the vector joining the centers of mass of 
the two molecules. The Hamiltonian of the system, 

H(r u r 2l R) = T(n) + T(r 2 ) + T{R) + V(r u r 2 , R), (1) 

is composed of a radial kinetic energy term T(i?) describ- 
ing the center-of-mass motion, two kinetic energy terms 



with 

Y\(r 1 ,h,R) = E (^i m A 1 A 2 mA 2 |Ai2mA 12 ) 

Y XimXl (h)Y X2mx2 (r 2 )Y^ mXi2 (R), (5) 

where A = A1A2A12 and m\ = mx^n^mx^. The 
notation (Xim\ 1 X 2 mx 2 \\i 2 m Xl2 ) represents a Clebsch- 
Gordan coefficient. Because the Hamiltonian defined in 
Eq. ([1]) is invariant under a rotation in space, the total an- 
gular momentum J and its projection J z on a space-fixed 
axis are conserved during the collision. In the following, 
we will discuss the close-coupling equations for both dis- 
tinguishable and indistinguishable molecule - molecule 
collisions. 

For collisions between two distinguishable molecules, 
the diabatic functions, 

<t> J v *{ ei (fi,f2,R) = X v jlj2 (ri,r2)(hr2R\jlJM), (6) 

serve as basis functions to expand the total wave func- 
tion, 

*(fi,f 2 ,^) = 

i E K™ S1 ( R )<t>t!? £I (ri,r%R), (7) 

v,j,l,J,M 

where v = viv 2 denotes the vibrational quantum num- 
bers and j = J1J2J12 denotes the rotational quantum 
numbers of the molecules. The function Xvj 1 j 2 { r ii r '2) = 
Xvtjt { r i)Xv 2 j2 ( r 2) is the product of the vibrational wave 
functions of the two diatomic fragments and the angular 
function, 

{hf 2 R\jUM) = 

E (h'm jl j 2 m n \j 12 m jl2 )(j 12 m n2 lmi\jM) 

rrij 1 ,rrij 2 ,rrij 12 ,mi 

Y Jim]1 {h)Y nm]2 {r 2 )Y lmi {R), (8) 
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denotes the rotational wave function of the two molecules 
in the total angular momentum representation. In this 
representation, the rotational angular momenta j\ and 
j2 are coupled to give the total rotational momentum 
ji2 which is subsequently coupled to the orbital angular 
momentum I to yield the total angular momentum J. 
Under spatial inversion /, the basis functions ((HJ) obey 
the eigenvalue equation, 



tiJMei _ i J Me i 

I<f> vjl = e I( /) vjl 



(9) 



with si = (-iyi+h+i = ±1. Substitution of Eq. ^ 
and Eq. (JT]) in the time-independent Schrodingcr equa- 
tion H*f> = E'fy leads to a set of close-coupling equations 
in the radial coordinate, R. The coupled equations be- 
come 



h 2 d 2 h 2 l(l + l) 



2/j, dR 2 2fiR 2 V1 

+ T,^mM r ) f ™i> 1 W = °> (10) 

v'j'l' 

where fi is the reduced mass of the molecule - molecule 
system. The total energy is E = e V j + E c where e v j = 
£ viji + £ii 2 j2 i s the rovibrational energy of the two sepa- 
rated molecular fragments and E c is the collision energy 
The interaction potential matrix is given by 



(11) 



The radial interaction potential term v *j'j'{R) is 

given by 



R vj!j 2 ,v 'j[ j' 2 (R) 



oo yoo 



Xvhj 2 (n,r 2 ) 

o Jo 

A\(r 1 ,r2,R)Xv'j' 1 j^(ri,r 2 )dr 1 dr 2 , (12) 

f ,/;A 



and the function fj^jip is given in terms of 3 — j, 6 — j, 
and 9 — j symbols by 

f-u'v = {^)- 3/ \-ir +h+ ^ +j [\3^f,i'M l/2 



3i fi AiW h 3 2 A 2 Ul V Aia 

o o o M o o oMoo o 



I V Ai 2 

j'\2 312 J 



3 12 32 3l 

312 h h } , (13) 
A12 A2 Ai 



with the notation 

] = (2xi + 1X2x2 + l)...(2ar„ + 1). (14) 

For two indistinguishable molecules, one has to sym- 
metrize the wave function with respect to the exchange- 
permutation symmetry of the molecules. Under molecule 
permutation P, P^f^ 1 = (-1)^+^+^+^^' , where 



v = V2V1 and j = 323x312- As a consequence, prop- 
erly symmetrized exchange-permutation invariant states 
must be employed for the basis set, 



' iJMei 



ep(-iy i+h+jl2+l cj)^ ei 



(15) 



with 



±1, A, 



[2(l + 5 VlV2 6 nn )]^ 2 , and 
= e P ^ JM£l£p 



The total wavefunction is expanded as 



(16) 



*(ri,f 2 ,-R) = 
1 

1; T, F - 



JMsiEf 



1£p (R)4> vj 

v,j,l,J,M 



JMeiei 



(fi,?t,R), (17) 



where the superscript over the sum symbol indicates a 
summation over states that have to follow a "well or- 
dered" classification as discussed in Ref . [iij [5(| [5l[ . We 
choose v\ > V2 or when V\ = v 2 , ji > 32- The time- 
independent Schrodingcr equation yields a set of close- 
coupling equations, 



H 2 d 2 h 2 l(l + \) 
^dR 2 + 2\iR 2 



E 



rpJMejEp 
vjl 



(R) 



l^Ktv'fv( R ) F v'7v iep ( R ) = 0, (18) 

v'j'l' 



with the interaction potential matrix, 



^vjl,v']'l'( R ) — ^vjin^v'^j'^yjiyyi'iR) 

+ E p(-l) j:+ii+j; ' +, '^, ( ,(fi)]. (19) 

Because of the indistinguishability of the molecules and 
the summation over the "well ordered" states, the size of 
the set of close-coupling equations (fP8)) is smaller than 
the set ([TO]) . 

The close-coupling equations (fTOf or (|18[> are solved 
for each value of R using the log -derivative matrix prop- 
agation method of Johnson 53 1 and Manolopoulos [54| . 
The log-derivative matrix is propagated to a matching 
distance where asymptotic boundary conditions are ap- 
plied to obtain the scattering matrix, S, for given values 
of J, si and sp. For distinguishable molecule collisions, 
the state-to-state cross section is given by 



E 

hsjigll'Jei 



(Ec) 



(2j 1 + l)(2j2 + l)k 2 



^ (2J+l)\S vjiyfl '-S^ v , j , l ,(E c )\ 2 , 



(20) 



where the wave vector k 2 = 2fiE c /h 2 . For indistinguish- 
able molecules, the state-to-state cross section is given by 
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a statistically weighted sum of the exchange-permutation 
symmetrized cross sections, 

where 

(2 J1 + l)(2j 2 + l)k 2 

(2J+l)IWfl'-^ P Wf. (22) 

For collisions between para-H2 molecules (nuclear spin 
I = 0), W + = 1,W~ = so that only the e P = +1 
exchange-permutation symmetry is needed in the close- 
coupling calculations, while for collisions between ortho- 
H 2 molecules (nuclear spin I = 1), W + = 2/3, W~ = 1/3 
and both exchange-permutation symmetries ep — ±1 are 
needed. The total inelastic cross section is the sum of all 
inelastic state-to-state cross sections. The state-to-state 
rate coefficients are obtained by integrating the corre- 
sponding cross sections over the Maxwell-Boltzmann dis- 
tribution of velocities, 



(T) = 



1 



8 



1/2 



{E c )e~ E ^^ E c dE c . (23) 



III. RESULTS AND DISCUSSION 

A. Computational details 

In this study we investigate the scattering of two 
identical para-H2 molecules using the BMKP PES [32| . 
For identical molecule collisions, the rovibrational quan- 
tum numbers obey the "well ordered states" classifica- 
tion [49|, [5(| [H| ■ For convenience we use the term "com- 
bined molecular state" (CMS), for a combination of two 



£.120' 




(0402) 
(0400) 
(0202) 
(0200) 
(0000) 



FIG. 2: The effective potentials, Vjr, as functions of R for 



the H4 system for J = 0, ej 



>eff> 

+l,e P = +1. 



rovibrational states of two H2 molecules. In this nota- 
tion, H2(ui, ji) + H2(v2, 32) collisions will be denoted as 
(W1J1W2J2) and its energy is given by e vj = e Vljl + e V2h . 
The CMS represents an unique quantum state of the di- 
atom - diatom system before or after the collision. To 
gain physical insights into the scattering process and also 
to roughly estimate the number of CMSs to be included 
in the scattering calculations, wc plot in Fig. [5] the effec- 
tive potentials which are composed of the diabatic ener- 
gies and the centrifugal terms, 



U 



JMejet 
vjl,vjl 



(R) 



h 2 i(i + i) 

2/j.R 2 ' 



(24) 



as functions of the radial coordinate, R, evaluated using 
the BMKP PES. At large separations, the energies of the 
different potential curves converge to that of the corre- 
sponding CMSs as indicated by the numbers (vijivzj?) 
in Fig'.CH The energy of the CMS (0000) is 6270.73 K and 
that of (0200) is 6780.35 K. For the vibrational^ excited 
case, (1000), the energy is 12257.57 K. The energy is rela- 
tive to the minimum of the H2 potential. Since the inter- 
action potential between two H 2 molecules is relatively 
weak the density of CMSs is rather sparse. The number 
of diabatic channels included in the calculation depends 
on the energy. For collision energies E c < 100 K, a cut-off 
energy is used to restrict the number of channels. This 
allows one to include all rotational and vibrational levels 
below the cut-off energy in the scattering calculations. 
We used cut-off energies of 20142.76 K (14000 cm" 1 ) for 
J = - 3, and 15826.45 K (11000 cm" 1 ) for J = 4 - 10 
for ultralow energy calculations. However, this approach 
may lead to a proliferation in the number of channels at 
higher energies as the cut-off energy needs to be higher. 
Thus only a restricted set of rovibrational levels are in- 
cluded in the calculations. In the present work we re- 
tained all CMSs with quantum numbers v = 0, 3 = — 8, 
and v = l,j = 0,2 for calculations at higher energies. 
The close-coupling equations (Tl8|) were propagated using 
the log-derivative matrix propagation method of John- 
son [53| and Manolopoulos [54| to an asymptotic match- 
ing distance of R = 53 do- For all collision energies, 
the inelastic cross section is converged at this value of 
the matching distance. However, to get converged elastic 
cross section for ultralow energies E c < 1 K, we propa- 
gate the elastic channel from R = 53 ao to a sufficiently 
large distance, at which the interaction potential U in 
Eq. ([2]) is less than 10~ 5 x E c . 



B. Vibrational relaxation in H2(u = l,j = 0) + 
T3.2(v = 0, 3 = 0) collisions 

In the upper panel of Fig. [3] we show the J-resolved par- 
tial cross sections and their sum for the elastic channel 
for collision energies ranging from 10~ 3 — 120 K. Though 
only the s-wave contributes at energies below 10 -2 K, the 
partial wave sum requires up to J = 10 at E c = 100 K to 
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yield converged cross sections. The state-to-state inelas- 
tic cross sections and their sum are shown in the lower 
panel of Fig. [3] for H 2 (u = 1, j = 0) + H 2 (u = 0, j = 0) 
collisions for E c = 10~ 3 — 120 K. The cross sections of the 
inelastic processes are much smaller than that of elastic 
scattering. Due to the small anisotropy of the interaction 
potential with respect to the stretching of the molecule, 
the interaction potential matrix elements responsible for 
vibrational relaxation channels are quite small compared 
to the isotropic part responsible for elastic scattering. 
Thus, the couplings between different CMSs are rather 
weak and state-to-state transitions are less probable than 
elastic scattering. The diabatic potential curves in Fig. [2] 
shows that the density of states is quite small and the 
energy levels are sparsely spaced leading to small in- 
clastic couplings. This is in contrast with the heavier 
alkali-metal systems with very high density of states, for 
which elastic scattering is less important than inelastic 
collisions at ultralow collision energies [l8||. In Fig. [3J 
the CMS (0800) is the most probable final state. This 
corresponds to a case where the vibrationally excited 
molecule H 2 (u = l,j = 0) relaxes to the rotationally 
excited = 0,j = 8) state upon collisions with the 

ground state molecule H 2 (v — 0, j = 0). The final state 
distribution is determined by a compromise between con- 
servation of the total internal energy and the rotational 
angular momenta of the colliding molecules. 

In Fig. [4] we present the state-to-state rovibrational 
populations of the two molecules in H 2 (w = 1, j = 0) + 
H 2 (v = 0,j = 0) collisions at an energy of E c = 10~ 3 K. 
Since rotational level resolved state-to-state cross sec- 
tions are more difficult to converge compared to initial 
state-selected cross sections, we also include results from 
our convergence studies. Figure [4] shows that the re- 
sults are robust and converged with respect to all nu- 
merical parameters employed in the calculations. The 
upper panel shows cross sections for different values of 
Ai = A 2 in Eq. It is seen that Ai = A 2 = 8 or 

higher is needed to accurately calculate cross sections 
for the dominant final CMS, (0800). The convergence 
of the results with respect to the cut-off energy, the size 
Ai? of the radial interval, the number of DVR points 
N ri = N r2 = N r ,N 9l = N 02 = N , and used for 
the vibrational and rotational wave functions are pre- 
sented in the lower panel. These results are obtained with 
Ai = A 2 = 8 for the angular expansion of the interaction 
potential. For the results reported in the final production 
calculations for the BMKP PES we used Ai = A 2 = 8, 
Ai? = 0.05 o , N g = 14, = 8, and iV r = 24. 

Figure 3] shows that vibrational relaxation of H 2 (u = 
1, j = 0) by collisions with H 2 (v = 0, j = 0) is driven by 
high order anisotropic terms of the BMKP PES. Pogreb- 
nya et al. [301 ] had previously investigated this issue and 
they had found that the high anisotropic terms of the 
BMKP PES lead to large values of the vibrational re- 
laxation rate coefficients compared to experimental re- 
sults. They found that a restricted version of this PES, 
referred to as the BMKPE PES which includes only 



H 2 (v=1,j=0) + H 2 (v=0,j=0) 




E C (K) 

FIG. 3: Upper panel: Partial-wave resolved and total elas- 
tic cross sections for H2(f = l,j — 0) + Bh(u = 0,,; = 0) 
collisions as a function of the collision energy. Lower panel: 
Inelastic state-to-state cross sections for the H2(« = 1, j = 0) 
+ H2(w = 0,j = 0) system as a function of the collision en- 
ergy. For clarity, only the final CMSs (OjOO) have been shown. 



Ai = A 2 = 2 components in Eq. ([4]) for the angular expan- 
sion of the interaction potential, yields better results in 
comparison with experimental data. In the upper panel 
of Fig. [5]we show the cross section for H 2 (t> = 1, j = 0) 
+ H 2 (u = 0,i = 0) -» H 2 (u = 0) + U 2 (v = 0) colli- 
sions for E c = 0.001 - 1 eV evaluated using the BMKPE 
(filled squares) and BMKP (filled circles) PESs. The 
corresponding results of Pogrcbnya et al. [1^ are also 
shown (open squares and circles) for comparison. While 
the results on the BMKPE PES are in overall good agree- 
ment, the results on the BMKP PES differ significantly at 
low collision energies (0.001 cV). The agreement improves 
with increase in collision energy and both cross sections 
show similar energy dependence for E c = 0.1 — 1 eV. In 
the lower panel of Fig. we compare state-to-state cross 
sections at E c = 0.01 eV for the BMKPE (white) and 
BMKP (black) PESs with the corresponding results of 
Pogrcbnya et al. (dashed or gray). Our results on both 
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H (v=1 ,j=0) + H,(v=0,j=0) ; E = 1 K 



H_(v=1 ,j=0) + H.(v=0,j=0) -> H.(v=0,all j) + H.(v=0,all j 
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FIG. 4: State-to-state cross sections for H2(v = l,j = 0) + 
Hh(« = 0,j — 0) collisions at an energy of 10 -3 K. The pa- 
rameters corresponding to the black distributions are adopted 
in the final production calculations. Upper panel: conver- 
gence with respect to the angular anisotropic terms of the 
PES. Lower panel: convergence with respect to the number 
of points for the radial and angular integrations and the cut- 
off energy for the basis set. 



PESs are in good agreement with those of Pogrebnya et 
al. for the CMSs (0000), (0200), (0202), (0400), (0402) 
but differ for the CMSs (0404), (0600), (0602), (0604), 
(0800). We believe that the differences are due to the 
CSA method employed in the work of Pogrebnya et al. 
This applies to the results in both panels of Fig. [5] The 
maximum value of the total angular momentum projec- 
tion quantum number, fl, on the body-fixed axis R used 
in the calculations of Pogrebnya et al. is |0| max = 4. 
In a body-fixed frame, the orbital angular momentum 
I is always perpendicular to the body-fixed axis R and 
the value mi of its projection on R is zero. This implies 
that the projection of J = ji2 + I on the body-fixed axis 
is O = rrij 12 . If |fi| max is restricted to 4, so does |m J12 |. 
This is equivalent to exclude all projection quantum num- 
bers |77ij 12 | > 4. Since no such restriction is imposed in 
our calculations, the differences between the exact and 
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FIG. 5: Upper panel: Vibrational relaxation cross sections 
as functions of the collision energy for the H2(« = l,j = 0) 
+ H2(u = 0, j = 0) system. Lower panel: State-to-state cross 
sections at a collision energy of E c — 0.01 eV. 



CSA methods for CMSs with quantum numbers ji2 > 4 
and \m jl2 \ > 4, such as (0404), (0600), (0602), (0604), 
(0800), are attributed to the restriction on fl imposed in 
the CSA calculations. The good agreement between the 
full quantum and CSA methods for the BMKPE PES 
is attributed to the negligible contributions of the final 
states (0404), (0600), (0602), (0604) and (0800). Since 
the restriction on O primarily impacts these final states, 
it does not significantly affect the cross sections computed 
using the BMKPE PES. 

The rate coefficients for the vibrational relaxation of 
H2(w = l,j = 0) in collision with H2(v = Q, j = 0) 
as functions of the temperature are reported in Fig. [S] 
for the BMKP PES (solid curve) and the BMKPE PES 
(dashed curve). The experimental results of Audibcrt 
et al. [46[ are also shown for comparison. For the rea- 
sons discussed above, the BMKP PES predicts rate co- 
efficients that are about two orders of magnitude larger 
than the experimental results for T < 200 K. The less 
anisotropic BMKPE PES yields results in better agree- 
ment with the experimental results. The agreement is 
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FIG. 6: (Color online) Rate coefficient for the vibrational re- 
laxation of H2(u = l,j = 0) in collisions withlL^ = 0,j = 0) 
as a function of the temperature. Results obtained using the 
BMKP and BMKPE PESs are compared with the experimen- 
tal results of Audibert et al. I4al. 



very good at T = 300 K. At lower temperatures, the 
calculations predict slightly larger values than the exper- 
imental results. For T > 300 K, the experimental result is 
higher than the theoretical results on the BMKPE PES. 



C. Rotational excitation and de-excitation in 

H 2 (v = 0, j = 0,2) + H 2 (w = 0,j = 0,2) collisions 

In Fig. [7| we show cross sections for H2(w = 0, j = 0) 
+ H 2 (w = 0,j = 0) (upper panel), H 2 (u = 0,j = 2) + 
H 2 (w = 0, j = 0) (middle panel) and H 2 (u = 0, j = 2) + 
H2(w = 0,j = 2) (lower panel) collisions on the BMKP 
PES. Previous results of Lee et al. [27[ on the BMKP and 
DJ PESs and Forrey [H] on the ZR PES based on the 
rigid rotor model are also shown for comparison. The re- 
sults on the BMKP PES in the upper panel show that the 
rigid rotor approximation is reliable in predicting elastic 
cross sections for collisions between two ground state H 2 
molecules. The DJ potential predicts similar results as 
the BMKP potential except for the slight shift in the 
location of the resonance at about E c = 2 K. The agree- 
ment between the DJ and BMKP results is excellent for 
collision energies above 4 K. Overall, the elastic cross 
section appears to be less sensitive to the choice of the 
PES and the dynamics approximation. In the middle 
panel of Fig. [7] we show the elastic and inelastic rota- 
tional de-excitation cross sections for H2(u = 0,j = 2) 
+ B.2(v = 0,j = 0) collisions. As in the case of colli- 
sions between ground state molecules, the elastic cross 
sections from rigid rotor calculations of Lee et al. [27| us- 
ing the DJ PES and the present full-dimensional results 
on the BMKP surface are in good agreement with each 
other. However, the inelastic cross sections show signif- 
icant differences. The inelastic rotational dc-excitation 
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FIG. 7: Elastic and inelastic cross sections as functions of 
the collision energy. Upper panel: Elastic cross sections for 
H 2 (« = 0,j = 0) + H 2 (u = 0,j = 0) collisions. Middle 
panel: Elastic and inelastic cross sections for H2(v — 0, j = 2) 
+ H2(u = 0,j = 0) collisions. Lower panel: Elastic and 
inelastic cross sections for H2(f = 0,j = 2) + H2(« = 0,j = 2) 
collisions. 



cross section on the DJ PES is about an order of mag- 
nitude larger than that derived from the BMKP PES. 
This confirms conclusions derived from previous calcu- 
lations [13, SI that the BMKP PES is not reliable in 
predicting accurate values of rotational excitation cross 
sections, especially for the Aj = 2 rotational transitions. 
In the lower panel we show elastic and inelastic rota- 
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tional dc-excitation cross sections for H2(u = 0, j = 2) 
+ H 2 (w = 0,j = 2) collisions. Interestingly, (0202) -> 
(0000) inelastic cross sections are comparable in magni- 
tude for the BMKP PES, the DJ PES, and also for the 
ZR PES. However, for the (0202) -> (0200) transition, 
the BMKP potential yields smaller cross sections than 
DJ and ZR potentials for collision energies below 2 K. 
The agreement improves for higher collision energies, es- 
pecially with that from the ZR potential. 

The elastic cross section of the H 2 (i' = 0, j = 0) + 
H 2 (i> = 0,j = 0) collision is shown in Fig. [8] for collision 
energies up to 0.03 cV. The theoretical results of Lee et 
al. [23 and experimental results of Bauer et al. [I?! are 
also included for comparison. The present results do not 
differ significantly from the results of Lee ct al. The full 
calculation predicts slightly larger values for the elastic 
cross section compared to the rigid rotor results. Overall, 
the DJ PES predicts results in better agreement with the 
experimental results. 
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FIG. 9: State-to-state cross sections for rotational excitation 
in H2(i; = 0,j = 0) + Hs(v = 0, j = 0) collisions as functions 
of the collision energy. The results of Otto et al. [4j§| are 
included for comparison. 



FIG. 8: (Color online) Elastic cross section for H2(i> = 0, j = 
0) + U2(v = 0,j = 0) collisions as a function of the incident 
collision energy. The experimental results of Bauer et al. [47| 
and theoretical results of Lee et al. [53] are also shown. 



A number of recent studies have reported rotational 
excitation cross sections in collisions between two rovi- 
brationally ground state H 2 molecules (27l . |42| . l43j . The 
studies of Gatti et al. [12] and Otto et al. [43[ employed 
the full-dimensional multiconfiguration time-dependent 
Hartree (MCTDH) approach while that of Lee et al. [27| 
employed the TIQM approach within the rigid rotor ap- 
proximation. Though the agreement between these cal- 
culations were generally good, important differences were 
found for some transitions, especially for (0000) — > (0202) 
and (0000) — > (0404) transitions. In Fig. [9] we show ro- 
tational excitation cross sections in Ha(« = 0,j = 0) + 
H 2 (w = 0, j = 0) collisions leading to final CMSs (0200), 
(0202), (0400), (0402), and (0404). To obtain converged 
cross sections, these calculations include contributions 
from J = - 100. Also included are TIQM rigid ro- 
tor results from the present work, and full-dimensional 



time-dependent calculations of Otto et al. (4j|. Figure 4 
of Ref. [431 presents similar comparison with the results 
of Lee et al. [2Z|, Lin et al. [IH, and Sultanov et al. [2gj |. 
Figure [9] shows that our results are generally in good 
agreement with those of Otto et al. for the (0000) — > 
(0200) and (0000) -> (0400) transitions. The agreement 
is better for collision energies lower than 0.5 eV. The rigid 
rotor calculations from the present study underestimate 
the cross sections. The largest discrepancy between the 
present work and the time-dependent calculations occur 
for the (0000) -> (0202) transition. Overall, for all of the 
transitions shown in Fig. [9] the agreement with the time- 
dependent results is better at low collision energies and 
it deteriorates with increasing collision energy. This also 
applies to the rigid rotor results compared to the full- 
dimensional calculations from the present study. Our 
full-dimensional calculations include v — 0, j — — 8 and 
v = 1, j = — 2 levels in the basis set. It is possible that 
a larger basis set with additional rotational levels would 
yield results in better agreement with the time-dependent 
results at higher energies. However, this does not explain 
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the discrepancy between the time-dependent and time- 
independent calculations for the (0000) — ► (0202) tran- 
sition. Calculations using significantly larger basis sets 
are computationally demanding and beyond the scope of 
this paper. 

H 2 (v=0,j=0) + H 2 (v=0,j=0) -> H 2 (v=0,j=2) + H 2 (v=0,j=0) 
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FIG. 10: Rate coefficient as a function of the temperature 
for U 2 (v = 0, j = 0) + H 2 (u = 0,j = 0) -» H 2 (« = 0, j = 2) 
+ H 2 (v = 0, j = 0) collisions. The TDWP results of Otto 
et af. [4^| and experimentaf resufts of Mate et af. [2f| are 
afso included for comparison. The TDWP results are aimost 
identicai to the present resufts for T > 10 K. 



In Fig. [TU]we show the temperature dependence of the 
rate coefficient for H2(v = 0, j = 0) + H2(v = 0, j = 0) — > 
H 2 (w = 0, j = 2) + H 2 (u = 0, j = 0) collisions along with 
the corresponding experimental results of Mate et al. [26| 
and the TDQW results of Otto et al. [H| . It has previ- 
ously been shown that for this transition the BMKP PES 
underestimates the rate coefficient by about an order of 
magnitude [13, EH . Our results confirm this finding. As 
Fig. llOl illustrates. our results are in very good agreement 
with those of Otto et al. [43[ although the correspond- 
ing cross sections given in Fig. [9] show some differences 
at collision energies above 0.5 eV. This is because in the 
temperature range of 10 — 300 K, the main contribution 
to the rate coefficients occurs from cross sections at en- 
ergies below 0.5 eV. 

IV. CONCLUSION 

In this paper we present explicit quantum scattering 
calculations of collisions between two para-H 2 molecules 



in v = and v = 1 vibrational levels from ultralow 
to thermal energies. A new code [48| based on a 
full-dimensional time-independent quantum scattering 
method has been employed for the calculations. Us- 
ing this code we have been able to test the reliability 
of the coupled-states approximation for accurately pre- 
dicting vibrational relaxation rate coefficients for H2 (v = 
l,j = 0) + H2(w = 0,j = 0) collisions. Comparison 
of our results with previous coupled-states calculations 
of Pogrebnya and Clary shows that the coupled-states 
approximation is not very reliable for this system in pre- 
dicting vibrational relaxation cross sections at low en- 
ergies. This is in part due to the high order angular 
anisotropic terms of the H4 BMKP PES. We obtained 
good agreement with the results of Pogrebnya and Clary 
for the BMKPE potential which excludes the high order 
anisotropic terms of the interaction potential. The re- 
laxation rate coefficient is found to be generally in good 
agreement with experimental results when the BMKPE 
PES is used. 

We explored pure rotational transitions in H2(u = 
0, j = 0, 2) + H 2 (v = 0, j = 0, 2) collisions at ultralow en- 
ergies and compared our results with those of Forrey and 
Lee et al., using different PESs. While the inelastic pro- 
cesses are sensitive to the PES employed, the elastic pro- 
cesses do not show significant differences. We also pre- 
sented pure rotational excitation of H2(v = 0, j = 0, 2) + 
H2(u = 0, j = 0, 2) and compared our results with those 
of Otto et al. The TIQM and TDWP methods give com- 
parable results for collision energies up to 0.5 eV. But 
some discrepancies exist for (0000) -► (0202) and (0000) 
— > (0404) rotational transitions, presumably due to the 
contributions from the v = 1 level. 

Future works will involve ortho-para and ortho-ortho 
collisions of hydrogen molecules as well as H2— HD, 
D2-D2 and HD-HD collisions. 
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